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O ' ABSTRACT 

o : 

■ We study both analytically and numerically hydrodynamical effects of two colliding shells, 
Oh. the simplified models of the internal shock in various relativistic outflows such as gamma-ray 

' bursts and blazars. We pay particular attention to three interesting cases: a pair of shells with 

the same rest mass density ("equal rest mass density"), a pair of shells with the same rest mass 

, ("equal mass"), and a pair of shells with the same bulk kinetic energy ("equal energtf) measured 

in the intersteller medium (ISM) frame. We find that the density profiles are significantly affected 

■ by the propagation of rarefaction waves. A split-feature appears at the contact discontinuity of 
l/"") , two shells for the "equal mass" case, while no significant split appears for the "equal energy" 

and "equal rest mass density" cases. The shell spreading with a few ten percent of the speed of 
^-j- . light is also shown as a notable aspect caused by rarefaction waves. The conversion efficiency 

of the bulk kinetic energy to internal one is numerically evaluated. The time evolutions of the 
efficiency show deviations from the widely-used incllastic two-point-mass-collision model. 

Subject headings: galaxies: jet — gamma rays: bursts — gamma rays: theory — radiation 
Q ' mechanisms: non-thermal — shock waves 

a: 

Ctf ■ 

1. INTRODUCTION 

•rH . 

The internal shock scenario proposed by Rees (1978) is one of the most promising models to explain 
the observational feature of relativistic outflows as in gamma-ray bursts, and blazars (e.g., Rees & Meszaros 
1992; Spada et al. 2001). In this scenario, the bulk kinetic energy of the outflowing plasma is converted 
into thermal energy and non-thermal particle energy by the shock dissipation and particle acceleration, 
respectively, and explain the large power of these objects. Based on this scenario, a lot of authors have 
attempted to link the observed temporal profiles to multiple internal interactions (e.g., Kobayashi, Piran, 
& Sari 1997 (hereafter KPS97); Panaitescu, Spada & Meszaros 1999; Tanihata et al. 2002; Nakar & Piran 
2002 (hereafter NP02)), looking for crucial hints on the central engine of these relativistic outflows. 

Most of the previous works focus on the comparison with the observed light curves and model 
predictions employing a simple inelastic collision of two point masses (KPS97) and little attention has 
been paid to hydrodynamical processes in the shell collision. However, it is obvious that, in the case of 
relativistic shocks, the time scales in which shock and rarefaction waves cross the shells are comparable to 
the dynamical time scale A /c, where A is the shell width measured in the comoving frame of the shell and 
c is the speed of light. Since the time scales of observations of these relativistic outflows (e.g., Takahashi 
et al. 2000 for blazar jet; Fishman & Meegan 1995 for GRBs) are much longer than the dynamical 
time scales, the light curves should contain the footprints of these hydrodynamical wave propagations. 
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Thus, it is very interesting to clarify the difference between the simple two-point-mass-collision (hereafter 
two-mass-collision) model and the hydrodynamical treatment. The recent study by Kobayashi & Sari 2001 
(hereafter KS01) reports that collided shells are reflected from each other by the thermal expansion. Since 
they perform a hydrodymamical simulation and show the reflection feature for a single case, the detail of 
propagations of rarefaction waves for various cases of collisions is not discussed. The aim of this paper 
is to clarify the hydrodynamical effects including the propagations of rarefaction wave. As the simplest 
case, we explore the hydrodynamics of two-shell-collisions in the internal shock model. Since we are mainly 
interested in the hydrodynamical processes themselves, it is beyond the scope of this paper to make a 
detailed comparison of the observed phenomena with the model results. 

We consider the time evolution of two colliding shells in relativistic hydrodynamics in §2. In §3, we 
discuss the application to GRBs and blazars. The summary and discussion are given in §4. 

2. HYDRODYNAMICS 

Here we consider the hydrodynamics of the two-shell-interactions. Our intention is to derive analytically 
various time scales for shocks and rarefaction waves crossing the shells. The fundamentals of relativistic 
shocks are given by Landau & Lifshitz (1959) and Blandford & McKee (1976). Our main assumptions are 

(1) we adopt a planar ID shock analysis and neglect radiative coolings for simplicity, (2) neglect the effect 
of magnetic fields, and (3)limit our attention to shells with relativistic speeds. We are currently planing 2D 
studies. The role of magnetic fields is still under debate. A multi-frequency analysis of TeV blazars shows 
that the energy density of magnetic field is smaller than that of non-thermal electrons (Kino, Takahara, & 
Kusunose 2002). As for (3), it is self-evident that the relativistic regime is most important, since emissions 
from GRB and blazars show a substantial Dopplcr boost. 

2.1. Shock Jump Condition 

In Fig. 1 we draw a schematic mass density profile during the shock propagation in the interactions 
of rapid and slow shells. Two shocks are formed: a reverse shock that goes into the rapid shell and a 
forward shock that propagates into the slow shell. There are four regions: (1) the unshocked slow shell, 

(2) the shocked slow shell, (3) the shocked rapid shell, and (4) the unshocked rapid shell. Thermodynamic 
quantities, shch as rest mass density p, pressure P, and internal energy density e are measured in the fluid 
rest frames. We use the terminology of regions i (i=l, 2, 3, and 4) and position of discontinuity j (j=FS, 
CD, and RS) where FS, CD, and RS stand for the forward shock front, contact discontinuity, and reverse 
shock front, respectively. The fluid velocity and Lorentz factor in the region i measured in the ISM rest 
frame are expressed as Vi(= Pic), and Ti, respectively. The relative velocity and Lorentz factor of the fluids 
in the regions i and j are denoted as Vij (= —Vji — ffijc — —/3jic) and IYj(= Ijj), respectively. Throughout 
this work, we use the assumption of I\ ^> 1. 

We first count the numbers of quantities and the shock jump conditions. Each region is specified 
by three physical quantities; rest-mass density pi, pressure P,, and velocity Vi measured in the ISM rest 
frame. Forward and reverse shock speeds measured in the frame of unshocked regions (i.e., regions 1 and 4, 
respectively) are two other quantities. In all, there are 3 x 4 + 2 = 14 physical quantities. The total number 
of the jump conditions at the forward shock (FS), reverse shock (RS), and contact discontinuity (CD) is 
3 + 3 + 2 = 8. Hence, given 3 + 3 = 6 upstream quantities for each shock, we can obtain the remaining 8 
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downstream quantities by using 8 jump conditions. 

Following Blandford & McKee (1976), we consider the limit of strong shock, P 2 /^2 3> Pi/ n i, and 
adopt the assumption that the upstream matter is cold. As an equation of state (EOS), we take 

Pi = (7i - l)(ei - Pi) (1) 
where ji is an adiabatic index. The jump conditions for the forward shock are written as follows: 

r2 (r 12 + i)[ 7 2(r 12 -i) + i] 2 

FS1 7 2(2 -72)(r 12 -l) + 2 ' 

p 2 72r i2 + 1 , s 

e2 = ri 2 p2 , — = — ; 7— , (2) 

Pi 72-I 

where Ti 2 = ^^(l — Pifo), and Tpsi is the Lorentz factor of forward shock measured in the rest frame 
of the unshocked slow shell. In the relativistic limit, the adiabatic index is 72 = 4/3. Using the same 
assumptions as in the forward shock, the jump conditions for the reverse shock are given by: 

2 _ (r 3 4 + i)[7 3 (r 3 4-i) + i] 2 

1 RS4 



7 3(2 - 73)(r 3 4-l) + 2 ' 
e 3 = 1 34P3 , — = — ; r— , (3) 

P4 73 -1 

where r 34 = r 3 r 4 (l — and Trs4 is the Lorentz factor of the reverse shock measured in the rest frame 

of the unshocked rapid shell. The equality of pressures and velocities across the contact discontinuity gives 

p 2 = p 3 , r 2 = r 3 . (4) 

Before a shock breakout, r 2 = T 3 = T^d is satisfied. It may be useful to introduce the ratio / = Pi/ p\, 
which can be obtained from P 2 = P 3 , as 

f = p± = (72ri2 + i)(r 12 -i) 

7 pi (7 3 r 3 4 + i)(r 34 - 1) ' u 

Throughout this paper, we set Pi = P4 = 0. Then, we take pi, p4, Ti, and T 4 as ramaining 4 upstream 
parameters. Then, with 8 shock conditions given above, we can obtain 8 downstream quantities p 2 e 2 , u 2 , 
«fs, P3, e 3 , v 3 , and v RS . 



2.2. Time Scales of Wave Propagations 

Here we evaluate seven time scales of relevance when the shock and rarefaction waves cross the colliding 
shells. They are useful to understand the hydrodynamical evolution of two-shell-collisions. We measure 
these time scales in the rest frame of CD (hereafter we call it the CD frame) because it facilitates the 
comparison with each other. In contrast, most of the previous papers used the ISM frame in measuring 
crossing time of shock (e.g., Sari & Piran 1995; Panaitescu et al. 1997). Here we need to introduce new 
physical paramaters, the shell widths measured in the ISM frame, A r , and A s , where subscripts r and s 
represent rapid and slow shells, respectively. In the ISM frame, the upstream parameters are as follows: 
Lorentz factors r r (— r 4 ), and r s (= Ti), rest mass densities p T , and p s . With 2 + 4 = 6 physical parameters, 
we can uniquely specify the initial condition. Note that the regions 1 and 4 disappear after FS and RS 
break out, respectively. 
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In the CD frame, we rewrite them as follows; A r , A s , T r (= T 4 = ^4), and T s (= T 1 = T12) during the 
shock propagation in the shells. After the shock breaks out of the shell, the velocity is not uniform and 
determined by the propagation of rarefaction wave. Note that once we choose the CD frame, T 4 and I\ are 
not independent of each other (see e.g., Eq. (8) below). 

The time in which FS crosses the slow shell, i FS , is given by 

A' 



t 



FS 



\Pl\ + \&s\ 



ri 2 r| sl (i-|/?;n/3 FS1 |)2 



1/2' 



(6) 



where we use Eq. (2) and /3 FS = 1 — l/r F 2 g . Thus, we can express i FS as a function of model parameters 
and T 1 which is also given implicitly from the model parameters. Similarly, the RS crossing-time in the 
rapid shell, t RS , is 



\&\ + Ks\ 



r; 2 r| S4 (i-|/3 4 ||/3R S 4l) 2 



1/2' 



(7) 



where we use Eq. (3) and /3^ s = 1 — l/Tjfg. It is important to note that T 1 and T 4 are not independent 
but related by the equation below, 



/(i - 73) + J.p(i - 7 3 ) 2 + 4/73 [(r; - lx^r; + 1) + /] 



2/73 



(8) 



It is expected that after FS has crossed the slow shell, a rarefaction wave (hereafter we call it FR) 
propagates into the shocked slow shell (e.g., Panaitescu et al. 1997). The sound speed is given by (e.g., 
Mihalas & Mihalas 1984) 



dP 

de 



ad e + F- 

Thus, the time at which FR reaches CD, £ FR _ CD , is given by 

A' 



t 



FR— CD 



FS 



*s,FS 
Cs2 



72-1 

7 2 r; + 1 



72(72 - i)(r; - 1) 

; ( 

72 r i -72 + I 



-1/2 



(9) 



(10) 



where A s FS is the width of the slow shell just after FS reaches the end of the shell. This is obtained by 
the mass conservation (e.g., Spada et al. 2001) where (72 — 1)/(72T 1 + 1) is the compression factor of the 
slow shell and T x is the factor from the Lorentz contraction. Just the same way, the corresponding time, 
i RR _ CD , at which the rarefaction wave (hereafter we call it RR) generated at the RS breakout reaches CD 
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is given by 



RR-CD 



*RS + 



-VRS 

Cs3 



%s + A r r z 



73 



7 3 r; + 1 



73(73 -i)(r;-i) t 

IS^'i - 73 + 1 



1-1/2 



(11) 



In the case of i RR _ CD > ^fr-cD' oiuv ^fr-cd ^ s an actual time and £ rr _cd ^ s a virtual time which does 
not exist in reality. The opposite case is also true. 

In the case of t RR _ CD < ipR-CDi we have the time at which two rarefaction waves collide each other, 

*RR-FR' as 



<-RR-FR 



t 



RR-CD 



*RR-CD + 



VRR-CD 



2c s2 

A s r i 72 — I / ^FR-CD ~ ^RR-CD 
2c s2 7 2 ri + 1 \ ipR-CD _ *FS 



(12) 



vhcrc A 



s, RR-CD 



is the width of the part of slow shell which FR has not passed through yet at t 



RR-CD' 



Since both RR and FR propagate at the speed c S 2 after £ rr _ce>, ^ ne above equation includes factor 2. Note 
that after the rarefaction wave crosses CD, the pressure gradient appears at CD. As a result, the CD begins 
to move from x = in the CD frame. Thus, Eqs. (12), (13), (14), and (15) are apploximated estimations. 
Similarly, in the case of i FR _ CD < ^ RR — cd 1 we have 



<-rr-fr 



t 



FR-CD 



+ 



A 



- t 



r, FR-CD 

2c s3 

A^r; 73-1 



RR-CD *FR-CD 



FR-CD 



2c s3 ^r; + 1 v *; 



RR-CD ^RS 



(13) 



The time scale in which RR catches up with the propagating forward shock (FS), i RR _FS' can be- 
estimated as 



t 



c s,2^ RR _CD 



RR-FS 



C S 2 - (3' FS c 

Similarly, the time-scale in which FR catches up with the reverse shock (RS) is approximated as 

Cr^FR-CD 



t 



FR-RS 



c s3 - /3 RS C 



(14) 



(15) 



2.3. Numerical Simulation 

We complcmentarily perform the special relativistic hydrodynamical simulations. The detail of the 
code is given in Mizuta et al. (2004). To sum up, the code is based on an approximate relativistic Riemann 
solver. The numerical flux is derived from Marquina's flux formula (Donat & Marquina 1996). This code is 
originally second order in space using the so-called MUSCL method. In this study, however, this is slightly 
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compromised for numerical stability. We assume plane symmetry and treat one dimensional motions of 
shells. In discussing the propagation of shock and rarefaction waves, we choose the CD frame. Given the 
ratio r r /r s in the ISM frame and the value of T s , we can determine the Lorentz transformation to the CD 
frame easily because the CD Lorentz factor T CD (= T 2 = T 3 ) measured in the ISM frame can be derived by 
solving Eq. (5). We should note that the number of free parameters is reduced from 6 to 5 because we have 
already fixed the frame. As for the EOS, we assume for simplicity that 73 = 72 = 4/3 for r 34 > 2, and 
73 =72 =5/3 otherwise. Although this simplification gives slightly inaccurate estimation on the speeds of 
wave propagations, there is little effect on our conclusions in this work. 

We start the calculation at t = when the collision of two shells has just begun. Thoughout this paper, 
we set A s /c = 1 and p r = 1 as units in numerical simulations. Initially, two shells have opposite velocities, 
namely, u r > and v s < 0. In the previous section, we did not impose any conditions for the plasma 
surrounding the two shells. We only assumed that the boundary of each shell will be kept intact during the 
passage of shocks and rarefaction waves. For our numerical runs, we put plasma of low rest mass density 
(1CP 4 <C 1, Ps/ Pr) outside of the two shells. They have the same velocity and pressure as the adjacent shell. 
At first, the boundary condition at the left boundary is a steady inflow of dilute plasma. When the reverse 
shock or the rarefaction wave reaches the rapid shell's boundary, the velocity of the dilute plasma is set to 
be zero instantaneously to reduce the effect of the interaction between the shock and the dilute plasma. At 
the same time, the left boundary condition is set to be a free outflow. The treatment of the right side dilute 
plasma and the boundary condition is the same as that of the left side. 

3. SHELL DYNAMICS AFTER COLLISION 
3.1. Shell Splitting 

3.1.1. General Consideration 

Here we classify the types of the mass density profile in the merged shell based on the order of the 
times obtained in the previous section. Table 1 gives the complete set of the possible orders. Although 
there are various cases in the orders, the density profile, in particular, the splitting feature is governed by 
the two criteria as follows. 

(I) When £r,r-cd < ^fr-cd f° r Pr > P$ or ^fr-cd < ^rr-cd f° r P$ > Pr is satisfied, the splitting occurs 
at the CD since the rarefaction wave going from the larger density region (region 2) into the smaller 
density region (region 3) makes a dip in the latter region. 

(II) When a pair of rarefaction waves propagating in the opposite directions collide with each other, the 
density begins to decrease at the collision point and the splitting feature emerges. Hence the existence 
of £rr-fr implies splitting feature. 

Based on these two criteria, the mass density profile is classified into four types and shown in Fig. 2. 
If both criteria are satisfied, then the mass density has triple peaks. We show the corresponding schematic 
picture of space time diagram in Fig. 3. If only one criteria is met, then the double-peaked profile is 
realized. When neither condition is satisfied, the single peak is obtained. 
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3.1.2. GRBs and Blazars 

Here we apply the above general consideration to the specific cases and examine which kind of rest 
mass density profile is realized in GRBs and blazars. We assume that the widths of two shells are same in 
the ISM frame which is written as A r /A s = 1 (sec, e.g., KS01). Wc consider following three cases since it 
seeems natural to suppose that ejected shells from the central engine have a correlation among them; (1) 
the energy of bulk motion of the rapid shell (E = Tmc 2 ) equals to that of the slow one in the ISM frame 
(we refer to it as "equal energy (or E)"), (2) the mass of rapid shell (to = pTA) equals to that of the slow 
one (hereafter we call it "equal mass (or m)"), and (3) the rest mass density of rapid shell equals to that of 
the slow one (hereafter we call it "equal rest mass density (or p)"), 

Px = Ps (equal p), 
p r A r r r = p s A s r s (equal to), 
p r A r r 2 c 2 = p s A s r 2 c 2 (equal £). ( 16 ) 

Note that in the case of A r = A s and I\ > T s , p s is always larger than p T . This leads to the absence of 
iRR-FS- For all cases, we have 3 + 1 + 1 = 5 parameters. We take the ratio r r /r s as the last parameter and 
vary its value. This completes the six model parameters. 

The various time scales for the "equal E" case in GRBs are shown as a function of r r /r s in Fig. 4. 
Here we set T s = 10 2 , A s = 10 10 cm, and p s = 10 -10 g cm -3 in the slow shell as an example. Slight jumps 
of time scales are seen at r r /r s ~ 5 in the Figure. They are caused by the abrupt change of adiabatic 
index between the non-relativistic regime and relativistic one. The softening of EOS in the relativistic 
regime leads to slower shock waves propagation in the CD frame. In the whole range, the criteria (I) and 
(II) given in the previous section are both satisfied. Therefore, the triple-peaked profile is expected (No. 
4 in Tabic 1) in principle. However, the criterion (I) is satisfied only marginally. As a result, two peaks 
are not remarkable. It is worthwhile to obtain order estimations of A r /A s and \(3 T /(3 S \ by using a simple 
apploximation of Tcd ~ T m (r m is the Lorentz factor of marged shell obtained by two-mass-collision model 
and it is given in the next subsection) in spite of some discrepancy with the exact solution of Eq. (5). We 
have 

a; _ A r r r rl ^ /r 4 \ 2 / r 2 + r| \ 
a; A s r s r; ~vrj Kn + n) 1 

As r r /r s increases, the ratios of shell widths and velocities in the CD frame go asymptotically to 

-3. (18) 

This equation explains well the fact that each time scale in Fig. 4 has a weak dependence on r r /r s . This 
is the reason why i RS and t FS are very close to each other. i RR „ CD is also close to ip R _cD simply because 
the sound speeds in the both shocked regions are about a few ten perccnts of the light speed and close 
to each other. The corresponding numerical results are shown in Figs. 5 and 6. In these calculations, we 
take the cases of r r /r s = 3 and r r /r s = 6, respectively. This implies a large density contrast of p s /pi = 9 
and p s / Pr = 36, respectively (see Table 2). The collision of the rarefaction waves occurs in the region with 
much lower density compared with region 2. As a result, the peak of the profile is smoothed out. For larger 
values of r r /r s , the density contrast between regions 2 and 3 becomes clearer. Hence we conclude that the 
"equal energy" case essencially evolves into single-peaked profiles. The space-time diagram obtained by the 



0T 



(r 2 + r 2 )(-r 2 + r 2 ) 



11)01 + r 2 ) 



(-r? 



(17) 



A, 
A' 
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numerical simulation for "equal £"' is shown in Figs. 5 and 6. From Fig. 5, we see that t FS ~ i RS ~ 3 as is 
shown in Eq. (18). In Fig. 6, these time scales become close to 3. 

The time scales for the "equal mass" case are shown in Fig. 7. Up to r r /r s ~ 20, the criteria (I) and 
(II) are both satisfied and the triple-peaked profile shows up (No. 6 in Table 1). As r r /T s increases, we 
obtain 



\ T4 



(19) 



and t RS and £rr_cd become larger compared with t FS and tpR-CD- For the numerical experiment, we 
select two cases which have r r /r s = 3 and 20. In each case, we clearly see the dip corresponding to the 
criterion (I) in Figs. 8 and 9. However, as in the "equal E" case, the collision of rarefaction waves occurs 
in the less dense rapid shell and the density peak tends to be smoothed out. Hence we conclude that the 
equal mass collision with a large value of r r /r s effectively evolves in to the double-peaked (D2) profile. In 
KS01, the authors found a shell-split feature in their numerical simulation (Fig. 2 in their paper) for the 
"equal m" case. It can be also explained as the D2 profile. In Fig. 9, we see that the rarefaction wave (FR) 
driven by the breakout of FS catches up with the shock wave (RS) from behind at t ~ 3, since the flow 
seen from RS is subsonic in the downstream of RS. The propagation speed of RS is modified by the merge 
with the rarefaction wave and is determined by the strengths of the shock wave and rarefaction wave. For 
the current case, the propagation speed of RS is almost unchanged up to its breakout at t £ 5.4. 

The time scales for the "equal p" case is shown in Fig. 11. The important point is that p s = p r along 
CD. Then the criterion (I) disappears. In the limit of large r r /T s , we have 



a. i m x 2 



a: 3 v ri 



i (20) 



and £rs and £rr— cd become larger compared with tps an d ^fr-cd- We see this in the numerical experiment 
with r r /r s = 6 in Fig. 11. As in Fig. 9, we see also in Fig. 11 that the FR catches up to the RS and 
merges during 3 <; t <^ 6.5. The last topic is the dependence of the above results on the hitherto fixed 
model parameters T s , A s , and p s . The results for different T s are completely the same as in Fig. 4. The 
position of each time scale is determined by A r /A s and P T //3 S and they are independent of A s itself. Hence 
we omit corresponding figures. For different p s , the result is also almost unchanged, since we treat the case 
of F r /r s > 3 for simplicity. Then the sound speed is about a few ten percent of the light speed and has 
a weak dependence on p s . When A s is increased, every time-scale increases linearly keeping the relative 
positions of the time scales. 

For a typical blazar, we set up three slow shell parameters as T s = 10, A s = 10 16 cm, and 
p s = 10~ 26 g cm~ 3 . We show the result in Fig. 12 for the "equal E" case as an example. The essential 
difference between GRBs and blazars is a typical shell width. Hence, as explained above, every time-scale 
becomes 10 6 times larger than that in Fig. 4 with the relative positions unchanged. 



3.1.3. Extra case 

Since we have not seen so far a clear case corresponding to the criterion (II), we have performed 
another specific case to show what happens for the collision of two rarefaction-rarefaction waves collision 
(Dl profile in Fig. 2). Dl profile is appears most clearly when the rapid and slow shells have similar mass 
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densities and shell widths in the CD frame. Hence we do not use the assumption of A r = A s here only. 
Instead, we employ the condition that p r = p s , A r = A s , and r r /T s = 20. In Fig. 13, Dl profile is indeed 
produced. A larger value of r r /T s produce a greater dip than that shown in Fig. 13. 



3.2. Shell Spreading 

In principle, we can obtain the speed of rarefaction wave using Riemann invariants (e.g., Zel'dovich & 
Raizcr, 1966). In the relativistic limit, it is known that the speed of the head of rarefaction wave is close 
to the speed of light (e.g., Anile 1989). As the EOS of the shocked region deviates from the relativistic 
one, the speed is reduced from the light speed and the intermediate regime should be treated by numerical 
calculations (e.g., Wen, Panaitescu, & Laguna 1997). It is worthwhile to note that from the values T r 
and T s in Table 2 we see that the EOS in the forward-shocked region is a non-relativistic one while the 
reverse-shocked region extends from the non-relativistic to the relativistic regime for the parameter ranges 
we adopted. 

Numerical results for shell spreading are shown in Figs. 5, 6, 8, 9, and 11, in which the width of the 
shell may be described as 

Kot(t) - K.FS + A r,RS + "FE(t' ~ 4s) + «Re(*' ~ 4s)> ( 21 ) 

where A tot is the total width of the shell in the CD frame. We should stress that although a lot of authors 
assume that the shell width is not changed after collisions (e.g., Spada et al. 2000; NP02) for simplicity, the 
reality is that the shells spread at vfe and wre in the forward and backward directions, respectively. 



3.3. Energy Conversion Efficiency 

The conversion efficiency of the bulk kinetic energy to the internal one is one of the most important 
issues to explore the nature of the central engine of relativistic outflows and many authors have studied it 
(e.g., Kumar 1999; Tanihata et al. 2002). 



3.3.1. Two-mass-collision model 

Let us briefly review the widely-used two-mass-collision model (e.g., Piran 1999). From the momentum 
and energy conservations, we have 

m r r r + m s r s = (m m + £ m /c 2 )r m , 
m r r r /3 r + m s r s /3 s = (m m + £ m /c 2 )T m f3 mi (22) 

where m m = m r + m B , £ m = £ T + £ s , and r m are the mass, the internal energy, and the Lorentz factor of 
the merged shell, respectively, and m r and m s are the rest mass of the rapid and slow shells, respectively. 
Then we obtain the efficiency e as 

e = 1 _ (m r + m s )r m r2 = r r m r r r + m s T s 
m r r r + m s T s ' m r s mr p s _|_ TOs p r 

It is a useful shortcut to approximate T m <~ T2 = T3 without solving Eq. (5). Using this shortcut, for 
the "equal mass density" case (m r /r r = m s /T s ), we have = (T 2 +T^)/2. For the "equal mass" case 
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(m r = m a ), we have = T T T S . For "equal energy" (m r r r = m s r s ), we have T 2 m = 2(r 2 r 2 )/(r 2 + r 2 ). 
Then, the efficiency e in each case is given by 



1/2 



4 



2+-^ + 



1 

72 



1 + 



(ft)'" 



-1/2 



T-l/2 



(equal p), 
(equal m), 
(equal £). 



(24) 



As the value of T r /r s gets larger for the "equal m" case, the efficiency becomes larger and approaches ~ 1. 
For "equal E" and "equal p" cases, it goes asymptotically to ~ 0.3. Thus we find that "equal m" is the 
most effective collision if we neglect the rarefaction waves. 



3.3.2. Shock model 

The consideration of shock dynamics provides us with the information of the assignment of the total 
dissipation energy E^ ss to the forward- and reverse-shocked regions Ep§ and -Ers- From the equations of 

mass continuity in the CD frame, we have A °T S = r } ^ — — and A "'" s = r . 4 ^j? — — Assuming a large value of 

A a 7 2 r 1 +i A r 7 3 r 4 +i 

r 4 /ri, then we have e 2 + Pi — e 3 + P 3 and the dissipated energies are mainly controlled by shell widths. 
They are written as 

A ' A ' 

EfS A s,FS ERS A r,RS ^5) 

-Bdiss ^ S ,FS + ^r,RS ^ diss ^s,FS + ^r,RS 

where Eps and Ers, internal energy of forward and reverse shocked regions just after the shock crosses 
each shell, respectively (e.g., Spada et al. 2000). However, the values E F s/E d [ BS and E R s/E d [ ss will begin 
to deviate from the the above approximation when the rarefaction waves start to propagate, since they 
reconvert the internal energy into the kinetic one. The result including the rarefaction waves is shown 
below. 



3.3.3. Numerical Results 

The estimation of the energy conversion efficiencies with shock and rarefaction waves taken into 
account are presented in Figs. 14, 15 and 16 based on the numerical simulations. From the analogy of the 
two-mass-collision model, we define the efficiency measured in the ISM frame as 

= 1 _ jT(t,x) dm(t,x) 
r r m r + T s m s 
_ fT(t,x) p(t,x)T(t,x)dx 

Pr A r r 2 + p s A s r 2 ' [D) 

where dm(t,x), p{t,x), T(t,x), and T{t,x)dx, are the rest mass element, the rest mass density, and the 
Lorentz factor, seen in the ISM frame, and the length of line element in the CD frame, respectively. We 
assume the origins of both frames coinside with each other at t = t = 0. The Lorentz transformation of 
T(t,x) is wtitten as T(t, x) — TqdT (t,x)[l — /3cd/3 {t,x)]. The hypersurface of a constant time in the 
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ISM frame does not coinside with that in CD frame. Hence, to evaluate spatial integrations at a certain 
time in the ISM frame, we must collect the values of physical quantities in the inegrand for different times 
in the CD frame. Unfortunately it is technically difficult task to perform this. Hence we report to an 
approximation that p(t, x) and T(t, x) are replaced by p(t ,x ) and T(t , x ), respectively. This is only valid 
near the original point of the CD frame and, elsewhere, mixes up those quantities at the different time slice. 
We believe, however, that this still gives the behavior of the efficiency and the essencial role of rarefaction 
waves. In these figures, we compare the numerical results with the prediction by the two-mass-collision 
approximation. We find that the two-mass-collision model well reproduces the hydrodynamical results just 
before rarefaction waves begin to propagate. After the shock waves break-out the shells, the conversion 
effeciency is reduced by several ten percent from the estimate of the two-mass-collision model after several 
dynamical times. It is noted again that no cooling effect is taken into account here. 

4. SUMMARY AND DISCUSSION 

In this paper we have illuminated the difference between the simple two-mass-collision model and 
the hydrodynamical one of the internal shock. We have studied ID hydrodynamical simulations of the 
two-shell-collisions in the CD frame taking the shock and rarefaction waves into account. Below we 
summarize our results and give some discussions. 

(1) By comparing the relevant time scales of shock and rarefaction waves, we have completely classified 
the evolutions of the two-shell-collisions using six physical parameters, that is, the widths, rest mass 
densities, and velocities of the two colliding shells. We find that rarefaction waves have a significant effect 
on the dynamics. In principle, the rest mass density profile can be evolved into single-, double-, and 
triple-peaked features. In the limit of r r ,r s » 1, the features are essentially characterized by only three 
parameters: the ratios of Lorentz factors, widths, and rest mass densities. The combination of the values of 
r r /r s and Pr//o s determines the relative orders of the time scales of various wave propagations, while the 
value of A r /A s controls normarizations of the time scales. 

(2) Bearing in mind the application to relativistic outflows such as GRBs and blazars, we specifically 
examine the cases of "equal p", "equal m", and "equal E" . For the "equal p" case, the profile is 
single-peaked. The rarefaction wave produced when the FS breaks out reaches CD and then catches up 
with RS. In the case of "equal m", the profile should in principle become triple-peaked according to our 
classification scheme. In practice, however, there is very little time for the FR-RR collision to make a 
clear dip, while there is a lot more time for the FR to create a dip for a fairly wide range of parameters. 
Therefore, the profile in this case is effectively double-peaked. For the "equal E" case, the profile is classified 
as triple-peak. However, again, there is very little time for the FR-RR collision to make a dip. A very large 
mass-density gradient between forward- and reverse-shocked regions makes the dip even less pronounced. 
Furthermore, there is again little time for FR to create a dip for a fairly wide range of parameters. Hence, 
we conclude that the profile for the "equal E" is effectively single-peakd. If the cooling time scale is 
sufficiently long in the shocked region, electromagnetic radiations will show these profiles. 

(3) For large r r /r s , we have shown that the spreading velocity of the shells after collision is close to 
the speed of light. Hence, the often used approximation of constant shell width after collision is not very 
good in treating multiple collisions. For examle, the authors in NP02 claim that the "equal energy" case is 
suggested for the shell Lorentz factors in GRBs, assuming that L > A, where L is the separation distance 
between two shells. If the interval of the first and the second collisions is long, however, the shell spreading 
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effect cannot be ignored and the case of L < A should be included in the analysis. Then the difference 
between the "equal to" and "equal E" cases might be wiped away. 

(4) As the shell spreads after collision, the internal energy is converted back to the bulk kinetic energy 
due to thermal expansion. We have numerically studied the timc-dependent energy conversion efficiency 
quantitatively. Since we have neglected cooling processes, the conversion effciency rises up to the order of 
unity. This should be corresponding to the event in the regime of "weak cooling" (Kino & Takahara 2004). 
If r r /r s ^> 1 and the time-interval between collisions is long, the conversion efficiency will be substantially 
deviated the estimate of the two-mass-collision model. 

We appreciate the insightful comments and suggestions of the referee. M.K. thank S. Kobayashi, K. 
Asano and F. Takahara for useful remarks and discussions. A.M. acknowledges support from Japan Sosiety 
for the Promotion of Science (JSPS). This work is also supported in part by the Grant-in- Aid Program for 
Scientific Research (14340066, 14740166 and 14079202) from the Ministry of Education, Science, Sports, 
and Culture of Japan. 
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Table 1. Various types of the evolution 



No. p timescale" profile 6 

1 Ps > Pi t RS < i RR _ CD < £rr_fs S 

2 ^RS < *RR-CD < *FS < ^RR-FR Dl 

3 ^RS < *FS < *RR-CD < ^RR-FR Dl 

4 ^RS < *FS < *FR-CD < *RR-FR T 

5 *FS < *FR-CD < *FR-RS 

6 i FS < t FR _CD < *RS < *RR-FR T 

7 i FS < t RS <FR-CD< *RR-FR T 
§ ^FS 5 ^RS < RR-CD < ^RR-FR Dl 

7 7 7 

9 Ps < Pr ^ t RS < i RR _cD < *RR-FS D2 

10 ^RS < *RR-CD < *FS < ^RR-FR T 

11 ^RS < *FS < *RR-CD < *RR-FR T 

12 ^RS < ^FS < *FR-CD < ^RR-FR Dl 

13 *FS < *FR-CD < *FR-RS ^ 

14 ^FS < *FR-CD < *RS < ^RR-FR Dl 

15 ^FS < *RS < *FR-CD < *RR-FR Dl 

16 ^FS < ^RS < ^RR-CD < ^RR-FR J_ 

17 Ps = Pr ^ t RS < t R R_ FS S 

1^ ^RS < *FS < ^RR-FR Dl 

19 ^ t FS < i FR _ RS S 

20 *FS < *RS < *RR-FR Dl 



a Notations are given in §2.2. 

b Representative profiles of S, Dl, D2, and T are schematically shown in Fig. 2. 
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Table 2. Parameter sets for numerical simulations of "equal E" , "equal m" and "equal p". 



No. a 


r r /r s 




r; 


r; 


K/K 


Ps/Pr 


4(equal E) 


3 


6.7 


1.34 


1.05 


2.3 


9 


4(equal E) c 


6 


7.0 


2.26 


1.06 


2.8 


36 


6 (equal m) 


3 


7.6 


1.25 


1.09 


2.6 


3 


5 (equal m) c 


20 


11.8 


4.29 


1.40 


6.4 


20 


19(equal p) 


6 


12.2 


1.43 


1.43 


6.0 


1 


-(equal p) c d 


20 


22.4 


2.35 


2.35 


1.0 


1 



a These numbers correspond to those shown in Table 1. 

b The values are obtained by solving Eq. (5). The slow shell Lorentz factor is fixed as T s = 5. 
c We use 73 = 4/3. 

d This corresponds to the case Dl profile in Fig. 2 without the assumption of A s = A r . 
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Fig. 1. — Sketch of a two-shcll-collision where a rapid shell catches up with a slower one at the CD frame. 
Forward and reverse shocks (FS and RS) propagate from the contact discontinuity (CD). We employ the 
conventional numbering for each region in the study of GRB (e.g., Piran 1999). 
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Fig. 2. — Schematic snapshot of the rest mass density of two shells. Here we assume p s > p r . The 
corresponding time scale relation in each case is shown in Table 1. Picture S shows the case where RR 
catches up with the propagating FS and has a single-peak. Picture Dl shows the case where both of FS and 
RS cross the shells and the rarefaction waves propagating from both sides dig a dip and give a double-peaked 
profile. Picture D2 shows the case where FR reaches CD and dig a dip. Hence the profile is double-peaked. 
Picture T shows the case of the combination of Dl and D2. Then the final profile is triple-peaked. 



- 18 - 




Fig. 3. — Schematic picture of space-time diagram of a two-shell-interaction in the CD frame. Here we 
assume p s > p T . A forward shock (FS) and a reverse shock (RS) run through the slow and rapid shells, 
respectively. After the shocks break out, the rarefaction waves propagate into the shells (FR and RR) and 
the shells spread (FE and RE). This corresponds to case 7 in Table 1. 
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Fig. 4. — r r /r s dependence of the various time scales for the "equal E"case. In the whole range, case 4 
(Triple) is realized. The slight jumps of the time scales at r r /r s <~ 5 correspond to the abrupt change of 
adiabatic index. The softening of EOS in the relativistic regime gives slower shock waves and the time scales 
get longer accordingly. 
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Fig. 5. — Left: Time evolution of the rest mass density profile in the CD frame in the case of "equal E" . 
The parameter is chosen so that r r /T s = 3 in the ISM frame. The parameters in the CD frame are shown 
in Table 2. Throughout our numerical simulations, we set A s /c = 1 and p s = 1 as units. Right: Space-time 
diagram of shock and rarefaction waves propagations. As is shown in the text, tus ~ *fs ~ 3. RE spreads 
at the speed ~ 0.9c while FE spreads at the speed ~ 0.6c. 




Fig. 6.— Left: Time evolution of the rest mass density profile in the CD frame for "equal E" . In the ISM 
frame, r r /T s = 6. The parameters are given in Table 2. Right: Space-time diagram of shock and rarefaction 
waves propagations. Both FS and RS propagate faster than those for r r /r s = 3 and ips ~ £rs ~ 3. RE 
spreads at the speed ~ c while FE spreads at the speed <~ 0.7c. 
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Fig. 7. — Time scales for "equal mass". When r r /r s is smaller than ~ 7, case 6 (Triple) is realized. For 
r r /r s larger than ~ 7, case 5 (D2) is realized. The slight jumps of the time scales at r r /r s ~ 7 correspond 
to the abrupt change of adiabatic index. The softening of EOS in the relativistic regime gives slower shock 
waves and the time scales get longer accordingly. 
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Fig. 8. — Left: Time evolution of the rest mass density profile in the CD frame for "equal to". In the 
ISM frame r r /r s = 3. The parameters are shown in Table 2. Criterion (II) is satisfied and the profile 
classified "D2" in Fig. 2 is seen in this numerical result. Space-time diagram of shock and rarefaction waves 
propagations. RE spreads at the speed ~ 0.8c while FE spreads at the speed ~ 0.7c. 




Fig. 9. — Left: Time evolution of the rest mass density profile in the CD frame for "equal to". In the ISM 
frame, r r /r s = 20. The parameters are shown in Table 2. The dilute rapid shell collides with the slow one 
and quickly spreads out. We also see that the dense slow shock is pushed forwards by the rapid shell. Right: 
Space-time diagram of shock and rarefaction waves propagations. RE spreads at the speed ~ c while FE 
spreads at the speed ~ 0.9c. 
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Fig. 10. — Time scales for "equal p" . Case 19 in Table 1 is realized in the whole range. The slight jumps of 
the time scales at r r /r s ~ 14, 15 correspond to the abrupt change of adiabatic index. The softening of EOS 
in the relativistic regime gives slower shock waves and the time scales get longer accordingly. 
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Fig. 11. — Left: Time evolution of the rest mass density profile in the CD frame for "equal p" . In the ISM 
frame, r r /r s = 6. The parameters are shown in Table 2. Right: Space-time diagram of shock and rarefaction 
waves propagations. RE spreads at the speed ~ c while FE spreads at the speed ~ 0.8c. 
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Fig. 12. — The same as Fig. 4 but for blazars. Compared with Fig. 4 the relative positions of various time 
scales are almost unchanged while the absolute value of each time-scale is reduced by a factor of 10 6 which 
is the ratio of the shell width of GRB to that of blazar. 
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Fig. 13. — Left: Time evolution of the rest mass density profile in the CD frame for "equal p" . In the ISM 
frame, r r /r s = 20. The parameters are shown in Table 2. By the rarefaction-rarefaction wave collision, the 
bump is generated at the center of the shell. This corresponds to the Dl profile in Fig. 2. Right: Space-time 
diagram of shock and rarefaction waves propagations. Both RE and FE spread at the speed ~ c. 
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Fig. 14. — Left: Time evolutions of the conversion efficiency defined by Eq. (26) for the "equal energy" case. 
During shock propagations, e approaches the two-mass-collision estimate given by Eq. (24) for each case. 
As the rarefaction waves begin to propagate, the efficiencies are reduced by them. This suppression is can 
be ascribed to the thermal expansion. 
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Fig. 15. — Time evolutions ol the conversion efficiency denned by Eq. (26) for the "equal m" case. During 
shock propagations, e approaches the two-point-mass estimate given by Eq. (24) for each case. 
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Fig. 16. — Time evolutions ol the conversion efficiency defined by Eq. (26) for the "equal p" case. During 
shock propagations, e approaches the two-point-mass estimate given by Eq. (24) for each case. 



